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ABSTRACT 

TIGER integral-field spectrography and Hubble Space Telescope WFPC2 imaging of 
the E3 galaxy NGC 2974 are used to derive the kinematics of the stellar and ionized 
gas components in its central 500 pc. We derive a numerical two-integral distribution 
function from a Multi-Gaussian Expansion (MGE) mass model using the Hunter & 
Qian formalism. The TIGER as well as published long-slit stellar kinematics, including 
higher order moments, are well fitted with this self-consistent model, requiring neither 
the addition of a significant mass contribution from a hidden disc structure, nor the 
presence of a central dark mass (at that spatial resolution). The data reveal the pres- 
ence of a striking, highly contrasted, two-arm gaseous spiral structure within a radius 
of ~ 200 pc, corresponding to a total mass of 6.8 x 10"^ solar masses of ionized gas. 
We use a deconvolved TIGER datacube to probe its kinematics at a resolution of about 
0'.'35 FWHM. Strong departures from circular motions are observed, as well as high 
velocity dispersion values on the inner side of the arms. We interpret the observed 
gas morphology and kinematics as the signature of streaming gas flows driven by a 
~ 540 pc diameter bar with Vtp = 700 ±100 km s~^ kpc~^. This hypothesis is strongly 
supported by the predictions of a density wave model. This model predicts that the 
bar should lie at about 35° from the line of nodes, and implies gas inflow towards 
the central ~ 50 pc. The quadrupole pertubation due to this bar is estimated to rep- 
resent less than 2% of the underlying gravitational potential (a maximum torque of 
about 10%), explaining the lack of a direct detection via broad-band photometry in 
the visible. Despite its weakness, the inner bar of NGC 2974 may be able to drive some 
gas within a 10 pc radius. We suggest that the presence of such inner bars might be 
more common among early- type disk galaxies than is generally thought, and that deep 
high-resolution emission-line imagery may be the best way to detect such structures. 

Key words: galaxies: individual: NGC 2974 - galaxies: nuclei - galaxies: interstellar 
matter, photometry, structure - galaxies: kinematics and dynamics 



1 INTRODUCTION 

It is now evident that early-type galaxies are far from be- 
ing the simple, (violently) relaxed, isothermal, purely stellar 
systems anticipated by the traditional picture developed by 
Hubble. They are often found to contain a complex inter- 
stellar medium (ISM), stellar discs and/or cusps and other 
discrete dynamical components (e.g. de Zeeuw et al. 2002). 
It now seems highly plausible that much of this complex- 
ity originates in mergers or other interactions between these 
galaxies and their environments. In addition, the statistical 
frequency of active galactic nuclei (hereafter AGNs) in the 
form of quasars and radio galaxies (which are generally as- 
sociated with elliptical galaxies) at medium to high redshift 
has led to the realization that the majority of nearby ellip- 



tical galaxies may well harbor a central massive black hole 
(e.g., Chokshi & Turner 1992). 

Indeed, recent stellar kinematic studies using Hubble 
Space Telescope (HST) spectroscopy have provided evidence 
for the presence of black holes in tens of early-type galaxies 
(Ferrarese k Merritt 2000; Gebhardt et al. 2000). They re- 
veal the presence of a correlation between black hole mass 
and bulge velocity dispersion, which suggests a causal con- 
nection between the formation history of the black hole and 
that of the host galaxy. 

One long-standing problem associated with the growth 
of black holes is the fueling mechanism. It has remained gen- 
erally unclear how to transfer significant amounts of mass 
from the central kiloparsec of galaxy into its inner regions 
(at a scale of 10 pc). While there are several known mecha- 
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nisms that can be invoked to initiate such gas inflow, most 
of them are commonly thought to be primarily relevant to 
spiral galaxies rather than early- type galaxies. 

One such mechanism is density waves, the presence of 
which can be inferred from structures like spirals or bars. 
These can develop naturally in galactic discs or be initiated 
by 'minor merger' of e.g. a massive spiral and a small satel- 
lite galaxy (e.g., Pfenniger & Norman 1990; Pfenniger 1991; 
Hernquist & Mihos 1995). In particular, inner 'secondary' 
bars in spiral galaxies (often referred to as 'bars within bars') 
have been identified as potentially important agents to fun- 
nel gas into the innermost tens of parsecs (e.g., Shlosman, 
Frank & Begelman 1989; Erwin & Sparke 1999; Emsellem et 
al. 2001). In this scenario, a large-scale bar transports mate- 
rial into a kiloparsec-scale disc. It may then be relayed by an 
inner bar that is able to transport the gas to within about 10 
pc of the galactic nucleus, approximately where the central 
supermassive black hole's potential can take over. Structures 
with a radial extent of 240 - 750 pc, interpreted as inner bars, 
are surprisingly common in early-type spirals (at least 25% 
of barred SO/Sa galaxies, and maybe as much as 40%, Erwin 
& Sparke 2002). Such small bars have sometimes been linked 
to the central activity, as in e.g. NGC 6946 where Elmegreen, 
Chromey & Santos (1998) mention the possible role of the 
210 pc bar- like structure and the central starburst. There is 
evidence that large-scale bars are responsible for an increase 
in the central gas concentration (Sakamoto et al. 1999). It 
is however important to emphasize that there is presently 
no strong direct evidence for the role of inner bars in gas 
fueling within the central ^ 50 pc. 

Martini & Pogge (1999) analysed visible and near- 
infrared HST images of 24 Seyfert 2 (spiral) galaxies and 
found that 20 of them exhibit spiral arms of dust within the 
inner few hundred parsecs. These inner spirals were thus sug- 
gested to play a crucial role in feeding gas into the central 
engines (see also Regan & Mulchaey 1999). The formation 
of such small gaseous spiral structures may involve acous- 
tic instabilities in non self- gravitating discs (Elmegreen et 
al. 1998 and references therein). However, the arm-interarm 
contrast of such spiral arms is typically less than 0.1 mag, 
and it is not yet clear how efficiently these spirals can drive 
gas towards the nucleus. A more important issue is that some 
of these spirals anyway require an external driver to exist, 
such as a tumbling potential. Emsellem & Ferruit (2000) re- 
ported on the discovery of such an inner gaseous spiral in 
M 104 (the Sombrero galaxy, a highly bulge-dominated Sa 
spiral galaxy) : HST imaging and integral-field spectroscopy 
hinted towards the presence of a (small) tumbling bar po- 
tential. 

Since most radio-loud AGNs are found in early-type 
(E and SO) galaxies, we started a HST program to study 
whether such inner spiral structures may also be present in 
early- type galaxies. We are targeting galaxies from the sam- 
ple of Goudfrooij et al. (1994a, 1994b), who performed a 
ground-based imaging and spectroscopic survey of ionized 
gas in a complete, magnitude- limited sample of 56 elliptical 
galaxies. We selected galaxies containing ionized gas fea- 
tures whose morphology and kinematics showed promising 
in terms of hosting contrasted nuclear gaseous spirals. In this 
paper we report on the results of HST imaging and integral- 
field spectroscopy of the E3 galaxy NGC 2974, where we 
discovered a strong inner trailing gas spiral. 



Table 1. Observational characteristics of the TIGER spectro- 
graphic exposures. The fields of view are only indicative as the 
merged TIGER exposures do not cover a rectangular area. 



Nov 1993 Run 


Blue 


Red 


Lens diameter 




0(^39 


Final spatial sampling 


0(^36 


0(^36 


Field of view 


8^4 X 7'/ 8 


7^4 X 7^8 


# of exposures 


3 


2 


Total exp. time 


135 mn 


105 mn 


# of merged spectra 


382 


429 


Spectral sampling 


1.5 A pixel" ^ 


1.5 A pixel"^ 


Spectral resolution 


4.0 A (FWHM) 


3.2 A (FWHM) 


Spectral domain 


5090 - 5580 A 


6550 - 6960 A 


Seeing (FWHM) 


1(^2 


0(^85 


Stellar template 


HR 1681 (KOIII) 





We first present the original datasets, which includes 
if5'T/WFPC2 imaging and TIGER integral-field spectroscopy 
(Sect. [21. We then proceed in Sect. with the two- 
dimensional maps resulting from a detailed analysis of the 
TIGER datacubes, and including a datacube deconvolution 
achieved with the use of the narrow-band HST imaging as 
a guide. The dynamical modeling conducted to fit the ob- 
served stellar and gaseous kinematics is presented in Sect.^] 
A brief discussion is provided in Sect.|5| followed by conclud- 
ing remarks in Sect. El 

In this paper, we will assume a distance of 21.5 Mpc 
(from surface brightness fiuctuations, worked out by Tonry 
et al. 2001), and use the variable d2i.5 = D [Mpc] /2 1.5 to 
normalize all quantities that are distance dependent. We 
favour this to a normalization like (e.g.) hrs (see Emsellem 
et al. 1996 for a rationale of this choice). This yields a scale 
of ^ 104 pc/arcsec. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 TIGER two-dimensional spectroscopy 

2.1.1 Data reduction 

We obtained integral field spectroscopy of NGC 2974 during 
a run in November 1993 at the Canada-France-Hawaii Tele- 
scope using the TIGER instrument. This spectrograph makes 
use of an array of about 500 micro- lenses whose diameter 
was set to 0^.^39 for these observations. Two configurations 
were used to cover the spectral domains around the MgA5172 
triplet (blue hereafter) and the Ha, [N ii]A6548,A6583 and 
[S ii]A6717,A6731 emission lines (red hereafter). The charac- 
teristics of both red and blue configurations are described 
in Table □ 

The reduction of these data was achieved with the 
TIGER software package written at the Lyon Observatory 
(Rousset 1992). This included bias subtraction, pixel-to- 
pixel fiat fielding, spectral extraction, wavelength calibra- 
tion, fiat fielding, removal of cosmic rays, differential atmo- 
spheric refraction correction and fiux calibration for all in- 
dividual exposures. The resulting data-cubes were then re- 
sampled to a square grid and merged. Details about these 
reduction steps can be found in Emsellem et al. (1996). In 
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the next paragraph, we only mention points specific to the 
Nov. 1993 TIGER run and the NGC 2974 data. 

A slightly incorrect position angle of the CCD with re- 
spect to the lens array led to some unwanted contaminations 
between adjacent spectra at their edges: this only affected 
the red configuration, for which spectra were safely trun- 
cated to a wavelength range of 6570-6960 A. The wavelength 
calibration was found to be accurate to within 0.03 pixel: 
this was made possible through the use of a Fabry-Perot 
etalon which provides regularly spaced Airy lines. The abso- 
lute flux calibration of the red spectra was performed using 
spectro-photometric standard stars. For the blue configura- 
tion we used an available flux- calibrated long-slit spectrum 
(Goudfrooij k Emsellem 1996). 

The stellar kinematics were derived from the blue 
spectra via a slightly modified version of Bender's (1990) 
Fourier Correlation Quotient (FCQ hereafter) method. Line 
strengths were measured using an empirical correction for 
the dispersion as in Emsellem et al. (1996). Since no spe- 
cific attempt was made to caUbrate our data onto the Lick 
system, we remind the reader that systematic offsets may 
exist between our line-strength values and other published 
ones. The emission-lines in the spectra were fitted by Gaus- 
sian profiles using the FIT/SPECTRA software written by 
Arlette Rousset (Lyon Observatory - Rousset 1992). 

2.1.2 Stellar template subtraction 

(a) The blue spectra: the [Ni]A5200 line 

As discussed in Goudfrooij & Emsellem (1996), when ionized 
gas is present the [N i] emission-line doublet at AA5198,5200 
(hereafter [Nl]A5200) can jeopardize the analysis of stellar 
dynamics and stellar populations of galaxies in the spectral 
domain around the Mgb feature . We therefore devised an it- 
erative procedure to eliminate the contribution of [N i]A5200 
to our blue spectra. This method will be described and il- 
lustrated in detail in a separate paper, so we only provide 
an outline here: 

(i) We first derive the stellar kinematics with FCQ and 
the original blue spectra as they are. 

(ii) We then use a library of pure absorption-line spec- 
tra, including template stars as well as galaxies devoid of 
emission lines, to fit each of our blue spectra, after proper 
broadening and redshifting (using the result of step (i)). A 
small spectral region around the expected location of the 
[Ni]A5200 line was masked out during this fitting process. 

(iii) We then subtract the resulting model absorption-line 
spectra and make a fit of the [N i]A5200 system on the resid- 
ual pure emission-line spectra. An example is presented in 
Fig.Q 

(iv) The fitted [N i] A5200 line systems are then subtracted 
from the original blue spectra, thus resulting in spectra from 
which the [Ni]A5200 contribution is eliminated. 

(v) Finally, we re-measure the velocities v and disper- 
sions a from the "cleaned" spectra: these values are the ones 
adopted for the remainder of the paper. 

(b) The red spectra 

Just as described above, we use a library of star and galaxy 
spectra to remove the stellar continuum in the red spec- 
tra. However, the limited wavelength coverage and the high 
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Figure 1. Example of the stellar template subtraction method 
for the blue spectra. The spectrum at the top shows the original 
spectrum (solid line) and best fit using a library of absorption-line 
template spectra (dotted line). At the bottom, the residual spec- 
trum (solid line) and the corresponding [Nl]A5200 line fit (dotted 
line) is presented. Flux unit is erg s"-*^ cm""^ arcsec"^ A"-*^. 
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Figure 2. Similar plots as in Fig. now for the red spectra. 
Two examples are shown, each time with the original spectrum 
(solid lines, top) and their respective fit (dotted lines, top), as 
well as the residual spectra (solid lines, bottom) and the Gaus- 
sian fit to the emission lines (dotted line, bottom). Flux unit is 
erg cm"-^ arcsec"-^ A"-*^. 

equivalent width of the emission lines in this spectral do- 
main makes this a somewhat difficult task. For that reason, 
we built an optimal template using the stellar kinematics 
derived from the blue spectra and only including the stellar 
templates which were used in the fit of the blue spectra: 
this prevents artificial variations of (e.g.) the width of the 
Ka absorption line. An example of such a fit and subtraction 
is shown in Fig. |21 The residual 'pure' emission-line spectra 
were then fitted using the fit/spectra software package 
(Rousset 1992) assuming Gaussian profiles for each individ- 
ual line. 



2.2 HST/WFPC2 imaging 

We used the Wide Field and Planetary Camera 2 (WFPC2) 
aboard HST to obtain high spatial resolution images of stars 
and ionized gas in the inner regions of NGC 2974 on April 16- 
17, 1997. Three exposures were taken with the Linear Ramp 
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Filter FR680 at a filter wheel rotation of +15°. This filter 
configuration results in a passband covering the Ha A6563 + 
[N II] AA6548,6583 emission hues at the redshift of NGC 2974 
when observed at the chosen position on the PC chip of 
WFPC2. Two exposures were also obtained with each of the 
filters F547M and F814W, whose passbands are free from 
(strong) emission lines at the redshift of NGC 2974. Total 
exposure times were 1400 s, 500 s, and 5100 s for the F547M, 
F814W, and FR680 filters, respectively. For each filter, one 
of the images was spatially offset by 0^.^5 from the others, 
which corresponds to an approximately integer pixel shift 
in both PC and WF CCDs. This was done to enable the 
identification of hot pixels during image combination. After 
standard pipeline processing and alignment of the images, 
we combined the images per filter using the STSDAS task 
CRREJ, which effectively removed both cosmic rays and hot 
pixels. Flatfielding of the FR680 images was done using a 
reference flatfield taken through the F658N filter. 

Since NGC 2974 contains prominent dust structures 
near its nucleus (see below and Fig. El, the removal of the 
continuum contribution from the FR680 image is more com- 
plicated than simply subtracting a F547M or F814W image. 
First of all, the F547M and F814W images were aligned with 
the FR680 image. The F547M and F814W images were then 
matched to the same PSF. As no bright stars are present 
on the WFPC2 frames, this was done by building PSFs for 
all three filters using tinytim (Krist 1992), and convolving 
each filtered image with the PSF of the other filtered image. 
A similar procedure was performed to match the PSF of 
the FR680 image to that of the F547M and F814W images. 
A synthetic image of the continuum at 6620 A was then 
constructed with the calcphot task in the SYNPHOT pack- 
age of STSDAS. We used the F547M and the ratio of the 
F547M and F814W images in conjunction with the assump- 
tions that (i) the unreddened spectral energy distribution 
of the stellar population matches that of the "E2" early- 
type galaxy spectrum of Bica (1988), and (ii) the dust fea- 
tures in NGC 2974 cause reddening that is described by 
the Galactic interstellar extinction curve of Rieke & Lebof- 
sky (1985), taking into account that half of the stellar light 
from NGC 2974 is in front of the nuclear dust features. This 
process will be described in more detail in Paper II of this 
series (Goudfrooij & Emsellem 2003, in preparation). This 
synthetic 6620 A continuum image was subtracted from the 
FR680 image, yielding the Ha+[Nii] emission- line image. A 
V—I colour map was also constructed, using the F547M and 
F814W images. Figure Elshows grey-scale representations of 
the F547M, F814W, V-I, and Ha+[Nii] images. 

These images will be discussed in detail in Goudfrooij 
& Emsellem (2003, in preparation) but will be used here for 
a datacube deconvolution (Sect. 13. 3 f) and in the modelling 
of the gaseous component (Sect.^^J- 



3 RESULTS FROM THE TIGER DATACUBES 

In the following Sections, we present the maps of the mor- 
phology and kinematics of the stellar and ionized gas compo- 
nents, as well as the Mgb, Fe5270 and Fe5335 line-strength 
maps. 

After barycentric correction, the systemic velocity of 
NGC 2974 was measured from the stellar absorption lines 



and found to be 1888 ±10 km s~^, significantly lower than 
the commonly used value of 1924 km s~^ (Davies et al. 
1987), but perfectly consistent with the value of 1890 it 
40 km s"^ quoted by Kim et al. (1988). All maps are oriented 
with the major axis along the horizontal axis. We remea- 
sured the position angle of the major axis in the outer part 
of the galaxy on direct TIGER images and found PA= 42.2° 
(the North axis will then be always 132.2 degrees from the 
vertical axis). Note that there is a slight isophote twist of 
about 4 degrees towards the centre, partly due to dust ex- 
tinction which is most prominent in the southeast side of 
the galaxy (cf. Fig.|^, thus defining its near side. 

In order to compare our results with previously pub- 
lished kinematics of NGC 2974 (using long-slit spectra), 
we simulated long-slit cuts through our data cubes follow- 
ing characteristics of the observations made by Bender et 
al. (1994) and Cinzano & van der Marel (1994; hereafter 
CvdM94) for the stellar component, and by Zeilinger et al. 
(1996) for the ionized gas component. Comparisons are plot- 
ted in Fig.2|and Fig.|51for the stellar and ionized gas compo- 
nents respectively, after subtracting the systemic velocity^. 
Both comparisons show excellent agreement. 

3.1 The stellar component 

The maps in this Section have all been reconstructed using 
the sky sampling of 0^.^36 for the pixel size: this way, one 
pixel in a map corresponds to one spatial element in the 
blue datacube. 

3.1.1 Stellar kinematics 

In Fig. |S] we present the stellar velocity V and dispersion a 
TIGER maps. The observed velocity field exhibits some de- 
partures from what is expected from axisymmetry. We mea- 
sure a tilt of the zero- velocity curve of 9° =1= 3°. This might 
partly be due to dust extinction perturbing the line-of-sight 
kinematics, although extinction is significant only on the 
South-East side. This should be confirmed by spectroscopy 
with higher signal-to-noise ratio. The central velocity gradi- 
ent is about 110 km s~^ arcsec"^ (or ^ 1094 d2i.5~^ km s~^ 
kpc~^). The velocity dispersion field is nearly flat in the cen- 
tral two arcseconds with a value around ~ 215 km s~^, and 
slowly decreases outwards. 

3.1.2 Line strengths 

We derived line strength maps from the blue TIGER spectra, 
namely Mgb, Fe5270, Fe5335 and Fe5406. Here we present 
the <Fe> map, (classically defined as (Fe5270 + Fe5335)/2) 
as well as the Mgb map before and after correction for the 
contribution of the [Ni]A5200 emission line (Figs. El and IHI, 
to illustrate its effect on the Mgb line strength (see Goud- 
frooij & Emsellem 1996). 

The highest Mgb value of 5.1 di 0.2 is reached about 
—0.9 arcsec from the centre along the major- axis, with its 
symmetric point having Mgb of 4.4il=0.2 (the central value is 
4.3 ± 0.2). There is a slight decrease of Mgb outwards, with 

^ We applied a barycentric correction to the Zeilinger et al. data 
sets taking into account the dates of their observations. 
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Figure 3. HST/WFPC2 images of the central 6x6 arcsec^ of NGC 2974. with a resolution of ~0('l (= 10.4 d2i.5 pc), except for the 
Hq;+[Nii] image in (d) which has resolution ^0^13. The orientation on the sky is indicated in panel (d). (a) The F547M (5454 A) image, 
and (b) the F814W (8269 A) image. The greyscales in both (a) and (b) use a logarithmic stretch. Darker shades represent brighter 
regions, (c) The V — I colour map, obtained from the ratio of (a) to (b) as explained in Sect. 12.21 The greyscales range linearly from of 
0.84 to 1.64 mag in V — I. Darker shades represent redder colours, (d) The Hq;+[Nii] image. The greyscales (in units of erg cm~^ s"-*^ 
arcsec"^) range logarithmically from 2.0 X 10~^^ to 4.0 x 10~^^, while the peak intensity is 9.6 X 10~^^. Darker shades represent brighter 
regions. 



values at the edge of the field being around 4. The <Fe> and 
Fe5406 maps are consistent with being flat. Higher signal- 
to-noise spectra would help confirming the absence of sig- 
nificant line strength gradients in the central arcseconds. 

3.2 The gaseous component 

We now present the emission-line maps derived from fits of 
the emission lines present in the red data cube. These maps 
were built using a pixel size of 0^.^18, half of the final sky 
sampling. 



3.2.1 Flux maps and spectra 



We detect Ha as well as the [N ii] and [S ii] doublets through- 
out the whole field of view. Ha is in fact observed on a much 
larger scale, as nicely emphasized with the maps produced 
by Plana et al. (1998). From our TIGER data, we cut out 
regions where the intensity of the [Nll]A6583 emission line 
was less than 5 times the 1 a noise level to avoid spurious 
features in the maps. This of course implies a lower flux 
threshold for the [S ii] lines which are significantly weaker 
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Figure 4. Comparison between the TIGER stellar kinematics 
(solid lines) for NGC 2974 with published data. Velocity V is 
shown in the left panels, and dispersion a in the right pan- 
els. Open squares are from CvdM94, and filled squares are 
from BSG94. Velocities are relative to the systemic velocity of 
NGC 2974. 



than the [N ii] lines. However, we kept the same field size for 
all lines to avoid any confusion. 

Spectral profiles of the emission lines show large line 
widths in the centre. This is true for both Ka and the 
forbidden lines ([Nil], [Sii]), so that it is not due to the 
presence of a Broad- Line Region of an active galactic nu- 
cleus (AGN). Spatial mapping of the ionized gas showing 
such broad lines shows that it is spatially unresolved in the 
TIGER data. We therefore modelled it with a fixed radial ve- 
locity Vwide and dispersion (Jwide- The values for ^^wide and 
o'wide were determined on the central spectrum, in which 
the wings have the best contrast and signal-to-noise ratio: 
we find Vwide — 1965 km s~"'^and (Jwide — 460 km s""*^. Er- 
rors bars on these values are difficult to derive since it should 
include both the formal error bar as well as the uncertainty 
due to the noise in the data. We estimated them using the 
25 spectra which are within 1^^ of the centre and were found 
to exhibit such a wide component. These central 25 spectra 
were fitted with an additional Gaussian component, and we 
then find rms uncertainties of 35 km s~^and 61 km s~^on 
the velocity and velocity dispersion respectively. All other 
spectra were fitted with a single Gaussian. 

The [Nii]A6583 emission- line flux map is presented in 
Fig.l^lwith spectra illustrating the variations of line profiles 



> 



200 


-200 

200 


-200 

200 


-200 



> 



-200 



J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l_ 

- VA □ 

135° 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i_ 

160° 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

□ 

a n n □ 
90° 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

1150 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

□ \ 
450 

J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

70° 


0° 

" ■ 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

:° "^^"^^ -- 

- - 

250 

" " 



-4 -2 2 4 -4 -2 2 4 
X [arcsec] x [arcsec] 

Figure 5. Comparison between the TIGER gas kinematics (solid 
line; Vg only) and the published long-slit kinematics (open 
squares) presented by Zeilinger et al. (1996). 



and hne ratios over the field. Other maps ([Sii] and Ha) 
are presented in Fig. 1101 They all exhibit a two-arm spiral 
structure that is consistent from one emission line to the 
other, with an additional bright central peak. The SW arm 
(negative abscissa) is brighter than the NE one. In the inner 
part of each arm, the emission-line profiles are complex. This 
is probably caused by the superposition of components at 
different velocities. 

Note that a simple a 'posteriori check of the valid- 
ity of the stellar continuum subtration method described 
in Sect. 12.1.21 can be done by reconstructing the two- 
dimensional [N i]A5200 intensity map and comparing it with 
the [Nii]A6583 map (Fig. (for which the surrounding 
continuum is essentially featureless) . The match is excellent 
(see Fig. EJ, particularly when considering the low signal- 
to- noise ratio of the residual [Ni]A5200 line spectra. 

The [Sii]A6717 / [Sii]A6731 ratio ranges from ^1 at 
the centre to ^1.4-1.5 in the arms, corresponding to elec- 
tron densities Ne of 500 cm~^ and less than 100 cm~^ (low 
density limit), respectively. The [Nii]A6583 / Ha ratio is 
roughly constant over the field of view, with values between 
2 and 2.2. This seems to contradict the measurements made 
by Zeilinger et al. (1996), who reported [Nii]A6583 / Ha 
values as high as 8 at the centre, a factor of two higher 
than ours at the centre. We suggest that this discrepancy is 
mainly due to the fact that Zeilinger et al. (1996) did not 
flux- calibrate their spectra (which were taken for the pur- 
pose of measuring gas kinematics) and did not correct them 
for the underlying (Ha) stellar absorption. 



® 2003 RAS, MNRAS 000.IT1IT71 



An inner bar in the early-type galaxy NGC 2974 ^ 





2 

X [arcsec] 



-2 2 

X [arcsec] 



-200 



V rkm/sl 



a rkm/sl 



260 



Figure 6. TIGER stellar kinematics of NGC 2974: surface bright- 
ness (top left), velocity (bottom left) and velocity dispersion (bot- 
tom right) maps. The step of the surface brightness contours is 
0.5 mag arcsec the velocity contour step being 20 km for 
both the velocity and dispersion. The central dots in each figure 
mark the centre of the isophotes. North has been orientated at 
132.2 degrees from the vertical axis as indicated. 
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Figure 7. Mg6 line strength maps before (left panel) and after 
(right panel) correction for the [NI]A5200 emission line. Surface 
brightness contours from the blue TIGER data (cf. Fig. [HJ are 
overplotted 
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Figure 8. <Fe>(left panel) and Fe5406 (right panel) line 
strength TIGER maps for NGC 2974. Contours as in Fig. 
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Figure 11. [Ni]A5200 fiux distribution map in NGC 2974. The 

step is 3 10"-*^^ erg s"-*^ cm~^ arcsec"^, and the brightest isocon- 
tour corresponds to 27 10"-*^^ erg s"-*^ cm"-^ arcsec"^. 



cleus is 1.9 x lO^^ ^^^^2 ^-^^^ a 4.7 x 10^^ dsi.s^ W 
contribution from the central unresolved peak. Assuming a 
uniform electron density A^e = 100 cm~^ and a pure hydro- 
gen gas, this corresponds to a mass of ionized gas in the arms 
of 6.8 X 10^ solar masses. For the central unresolved peak 
and assuming this time an electron density of 500 cm~^, we 
obtain a mass of ionized gas of 4.4 x 10^ solar masses. Note 
that the quoted ionized gas masses are inversely proportional 
to the value of Ne. 



The flux distribution of the [S ii] lines is less cen- 
trally peaked than that of the [N ii] lines, with the 
[Sii]AA6717,6731 / [Nii]A6583 ratio ranging from :^0.4 at 
the location of the nucleus to ^0.65 in the arms. This is 
usually interpreted as a decrease of the [S ii] emission in the 
dense nuclear regions, due to the lower critical electronic 
densities of the [S ii] lines with respect to the [N ii] line dou- 



blet 



1.9 X 10^ 1.2 X 10^ and 8 x 10^ cm"'' for the 



[Sii]A6717, [Sii]A6731 and [N ii]AA6548,6583 lines, respec- 
tively) . 

The Ha flux integrated over a 3^^ radius from the nu- 



3.2.2 Gas kinematics 

The emission-line velocity and velocity dispersion maps are 
shown in Fig. 1121 The velocity field exhibits strong distur- 
bances following the spiral-like distribution of the gas. There 
is a strong hint for gas streaming on the inner side of the 
south-western arm, where the emission lines exhibit a com- 
plex structure probably resulting from the superposition of 
several kinematical components (Fig. This widening of 
the lines (dgas > 250 km s~^) is also observed on the inner 
side of the north-eastern arm, as emphasized in the disper- 
sion map (Fig. 
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Figure 9. Flux map of the [NII]A6583 emission line derived from the red TIGER datacube. Spectra (solid lines) and their corresponding 
multi-Gaussian fits (dashed lines) are shown in various places over the field. Note that all fiux scales are the same, except for the central 
spectrum shown in the central top panel. Flux unit is 10"-*^^ erg s"-*^ cm~^ arcsec"^ for the spectra. 



3.3 Datacube deconvolution 

Using the high-resolution if5'T/WFPC2 narrow-band 
Ha+[Nii] image of NGC 2974 shown in Fig. El we applied 
the method pioneered by Ferruit et al. (1999) to deconvolve 
our merged TIGER datacube (see also Emsellem & Ferruit 
2000). Two methods were tested: a pure Lucy deconvolu- 
tion, and a "weakly guided" Lucy deconvolution in which 
the narrow-band if5'T/WFPC2 image is used to constrain 
the integrated flux in the red datacube. In both cases, we 
performed 300 iterations after which the gain in resolution 
was found not to compensate the increase in noise level. Dif- 
ferences between the results of the two deconvolutions were 
not significant, except for a slight reduction of the high- 
frequency features in the weakly guided case. We therefore 
favour the latter, the analysis of which will be presented 
here. 

The final resolution of the deconvolved datacube, eval- 
uated from the central peak (unresolved in the WFPC2 
[Nii]+Ha image), is 0'.'35 FWHM. Maps of the deconvolved 
datacube, namely the [N ii] and Ha flux distribution, the gas 
velocity and dispersion maps, are shown in Fig. 1131 

The two spiral arms are nicely revealed, as well as 
the central concentration which is now highly contrasted. 
The velocity field shows strong streaming motions along the 
arms, with peak velocities of —211 and 264 km s~^. More 
interestingly, the dispersion is very high on the inner side of 



both arms, confirming the picture seen before deconvolution 
(see Fig. IT^ . reaching values greater than 300 km s~^. The 
gas kinematics will be further discussed in Sect. 14.31 in view 
of a density wave model. 



4 DYNAMICAL MODELS 

In this Section, we present new dynamical models for 
the stellar and gas components in the central region of 
NGC 2974, using our new two-dimensional data as well as 
published kinematics. 



4.1 The MGE luminosity model 

The first step in this modelling process is to obtain a re- 
alistic three-dimensional representation of the luminosity 
and mass distribution of the galaxy from the very centre 
to the outer region. We used the Mult i- Gaussian Expan- 
sion method (Monnet, Bacon & Emsellem 1992; Emsellem 
et al. 1994) to build a photometric modeP for the decon- 
volved surface brightness of NGC 2974. the combination of 
a ground-based / band image taken at the 1.0-m JKT (data 



^ See Cappellari 2002 for an efficient MGE package developed 
within the IDL environment. 
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Figure 10. Flux maps of [Nll]A6583 (top left), Ha (top right), [Sll]A6717 (bottom left) and [Sll]A6731 (bottom right). The lowest 
isocontours corresponds to lO"-*^^ erg s"-*^ cm"-^ arcsec"^, and the isocontour step is lO"-*^^ erg s"-*^ cm~^ arcsec"^ for [Nil] and Ha, 
and 5 10~^^ erg s"-*^ cm~^ arcsec"^ for both [Sll] maps. 
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Figure 12. Gas velocity (left) and velocity dispersion (right) maps for the emission lines in the red datacube. The isovelocity con- 
tours have a step of 50 km . The isocontours superimposed on the dispersion map correspond to the [Nll]A6583 flux map (step of 
lO"-*^^ erg s"-*^ cm~^ arcsec"^). 
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Figure 13. Maps built from the deconvolved TIGER datacube: [Nll]A6583 and Ha flux distribution (top left and right panels respec- 
tively), gas velocity (bottom left) and velocity dispersion (bottom right) maps. The step is 10"-*^^ erg s"-*^ cm~^ arcsec"^, and the 
lower level contours correspond to 60 and 30 erg s"-*^ cm~^ arcsec^-^for the [Nil] and Ho; lines respectively. The isovelocity con- 
tours have a step of 30 km s"-*^. The isocontours superimposed on the dispersion map correspond to the [Nll]A6583 flux map (step of 
lO"-*^^ erg s"-*^ cm~^ arcsec"^). 



taken from Goudfrooij et al. 1994a) and the if5'T/WFPC2 
F814W image provided the required wide field of view and 
high spatial resolution, respectively. The result of the fit, 
taking into account the respective Point Spread Functions 
of the input images, is shown in Fig. 1141 Note that the re- 
gions affected by dust extinction were masked during the 
fitting procedure, although we checked that this did not af- 
fect the result of the fit in any significant way. The spatial 
luminosity density is obtained by assuming that each in- 
dividual Gaussian component is uniquely deprojected as a 
three-dimensional axisymmetric Gaussian (see Table 

The total luminosity of the model is 2.9 10^^^21.5^ LqIu 
the / band. The most fiattened gaussian component in the 
deprojected model is Gaussian #7, with an axis ratio of 
0.344, which contributes to about 19% of the total mass in 
the model. 



# 


/(L0.pc ^.arcsec ^) 


an 




1 


807834.22 


0.048 


0.782 


2 


355349.63 


0.105 


0.608 


3 


62166.935 


0.313 


0.539 


4 


16377.007 


0.605 


0.649 


5 


3738.1343 


1.326 


0.598 


6 


627.7568 


3.735 


0.501 


7 


104.9081 


9.667 


0.344 


8 


9.3040 


21.829 


0.461 


9 


0.8656 


55.620 


0.466 



Table 2. Deprojected parameters for the MGE luminosity model 
of NGC 2974. From left to right: component number, maxi- 
mum deprojected intensity, width as measured by the Gaussian 
sigma, and axis ratio. The inclination angle was fixed to 60° (see 
Sect.lOI 



4.2 Two-integral models for the stellar component 

As detailed in Emsellem, Dejonghe, & Bacon (1999, here- 
after EDB99), we constrained the inclination and mass-to- 
light ratio by using simple dynamical Jeans modelling. A 
best-fit model was obtained with M/Li — 4.6 d2i.b~^ and 
an inclination oi i — 60°, consistent with the inclination 
found by Plana et al. (1998). However, all models with val- 
ues 58 < i < 65 were found to be consistent with the ob- 
served kinematics. No central dark mass was required to fit 
the central velocity gradient and velocity dispersion value at 



the resolution of our data (see Fig. I15f) . The observed kine- 
matics imply an upper limit of 2.5 lO^c/21.5 to the central 
dark mass. This is consistent with the mass inferred using 
the Mbh/c" correlation (Ferrarese & Merritt 2000, Gebhardt 
et al. 2000). 

We then derived the even part of the two-integral distri- 
bution function fe{E,Lz) corresponding to the MGE mass 
model via the Hunter & Qian (1993) routine (details of 
this procedure can be found in EDB99). The odd part fo 
was parametrized with respect to the even part using the 
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Figure 14. Comparison between our MGE model and photometric data in the / band. The left-hand panels depict the JKT data, and the 
right-hand panels depict the WFPC2/F814W data. Top panels: I-band isophotes (thin contours; step is 0.5 mag) and the corresponding 
MGE model isophotes (thick contours). Two bottom panels: major- and minor-axis profiles of the JKT and WFPC2/F814W images 
(thin lines) and of the MGE model (thick lines). 
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Figure 15. Comparison between the stellar kinematics of 
CvdM94 along the major (top) and minor (bottom) axis and 
isotropic Jeans models with different black hole masses (solid, 
dotted and dashed lines for 0, 3 and 5 10^ M©, respectively). 



function ha{r]) [where r] = J/Jmax{E)]\ E.g., fo{E,J) = 
haivi) fe{E, J) (see Eq. 1 in EDB99, and references therein). 
The free parameter a was taken to be the same for all com- 
ponents. 

Line Of Sight Velocity Distributions (LOSVDs) corre- 
sponding to the computed distribution function were then 
derived on a square grid, and convolved by the appropriate 
kernel to include the effect of seeing and pixel integration. 
These LOSVDs were parametrized via Gauss-Hermite mo- 
ments as well as true velocity moments, and compared with 
the observed kinematics. 



4.2.1 The best- fit model 

Our best-fit model has M/Lj = 4.63 (^21.5"^ with a = 2.5 at 
the centre, and a = 5 for ^/^"^"^ < 0.38 corresponding to 
a circular orbital radius of Rc > The agreement with the 
CvdM94 long-slit data is excellent; a comparison is shown in 
Fig. Although the model fits the TIGER maps relatively 
well, it can obviously not reproduce the slight tilt of the 
zero velocity curve (Fig.ini. Our value for M/Li is larger 
than the one found by CvdM94 (who found M/Lr '-^ 5.19 or 
M/Li ^ 3.92, normalized to our assumed distance). How- 
ever the total mass of the galaxy, as measured on the MGE 
model, is 1.34 10^^^21.5^ consistent with the value de- 

rived by CvdM94 (again for the present assumed distance of 
21.5 Mpc). We also find an excellent agreement between our 
value for the predicted circular velocity at 2\ 219 km s~^, 
and the one derived by CvdM94, 217 km s"^ This means 
that this mass to light ratio discrepancy comes mostly from a 
difference of zero points in the photometry of about 0.2 mag- 
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Figure 16. Comparison between our f(E, Lz) model (solid lines) 
and the stellar kinematics published by CvdM94 (same symbols 
as in CvdM94). 



nitude. Our estimate is however significantly higher than the 
one obtained by Pizzella et al. (1997) for a triaxial geometry. 

Fig. 1181 provides a representation of the total best-fit 
distribution function (even plus odd parts) versus the two 
integrals of motions. This plot can be compared to the one 
derived for NGC 3115 (EDB99). In contrast with NGC 2974, 
NGC 3115 posseses rather bright and very thin disc struc- 
tures, which can be easily identified in the even part of the 
distribution function as strong peaks near the expected lo- 
cation for the circular orbits (J/J"^""^ = !)• 



4.3 Gaseous Density wave modelling 

As mentioned in Sect. 13.2.11 the observed two-dimensional 
morphology and velocity field of the ionized gas in the cen- 
tral few arcsec show the presence of a highly contrasted, 
two-arm spiral structure with strong streaming motions. At 
a radius of 2'\ corresponding to about 220 pc at the distance 
of NGC 2974, and for a velocity of about 200 km s"\ the 
dynamical timescale is only about 10^ years. The contrast 
and large opening of the spiral arms are therefore suggesting 
that a quasi- stationary process is at play here. In what fol- 
lows, we thus present our attempt to fit the gas morphology 
and kinematics with a density wave model. 



4.3.1 The geometry of the spiral arms 

We first need to constrain the parameters of the spiral arms, 
namely its geometry and amplitude. This was done with 
the help of the Ha+[Nii] WFPC2 image corrected for ex- 
tinction (Goudfrooij & Emsellem 2003, in preparation). The 
narrow-band image was deprojected assuming an inclination 
angle i — 60° as obtained from the stellar dynamical model 
(Sect.2I3>- We then estimated the pitch angle as a function 



© 2003 RAS, MNRAS QQQ.IT1IT71 



An inner bar in the early-type galaxy NGC 2974 13 




-4-2 2 
X [arcsec] 



-4-2 2 
X [arcsec] 



Figure 17. Comparison between our f(E,Lz) model and the 
TIGER stellar kinematics. Top: observed (left) and modeled veloc- 
ity fields. Bottom: major (solid lines) and minor (dashed lines) 
axis velocity (left) and dispersion (right) profiles. The thick lines 
correspond to the model, the thin lines to the observed kinemat- 
ics. The minor-axis dispersion profile has been shifted for clarity. 
Note the slight minor-axis velocity gradient, which is not repro- 
duced by the axisymmetric model. 
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Figure 19. Pitch angle of the ionized gas spiral between a ra- 
dius of 0^8 and 2^^, as measured from the deprojected Hq;+[Nii] 
WFPC2 image. 



source of emission dominates (see Fig. I21|) . We can follow 
the spiral up to a deprojected radius of about 0V8 (i.e. ap- 
parent radius of 0^4, as the connection of the spiral with the 
central source unfortunately occurs close to the minor-axis 
of the galaxy). Starting at a radius of lV5, an abrupt change 
in the pitch angle occurs: it decreases from 40° to almost 0° 
in 0V5. 




Figure 18. f(E, Lz) model represented as an illuminated surface 
of log ((if) versus the normalized energy E / E^^^ and angular 
momentum J j J^""^ . 



of radius. The gaseous arm on the NE side of the centre is 
better defined and less patchy, probably due to the lower ex- 
tinction (the near side of the galaxy being the SE side). We 
therefore decided to preferentially use the NE arm to mea- 
sure the pitch angle. The average pitch angle is around 35° 
(see Fig. with values between 20° and 50°. It seems to 
increase towards the centre, although the spiral is not easily 
traceable close to the nucleus where the unresolved central 



Hints for the presence of a tumbling potential 

The total mass of ionized gas within the central 3^^ is ^ 
6.8 10^ M0, to be compared with a stellar mass of 

0.8 10^° d2i.5^ Mo derived from the MGE model (Sect.g^. 
We can safely say that the gas component is therefore not 
self- gravitating in this region. In the limit of tightly wound 
spirals, we could use the dispersion relation for non self- 
gravitating m — 2 spirals to derive the pitch angle (assuming 
a value for the sound speed Vs^ see e.g. Englmaier & Shlos- 
man 2000). Such density waves are supported inside their 
own Inner Lindblad Resonance (ILR hereafter), or outside 
the Outer Lindblad Resonance (Maciejewski et al. 2002). 
However, the large pitch angle observed in the case of the 
inner gaseous spiral arms of NGC 2974 prevents us to use 
such an approximation. 

More importantly, we can also exclude the role of 
self- gravity as the source for the observed m = 2 mode. 
We therefore need to search for another driver to explain 
the highly contrasted set of arms in the inner 300 pc of 
NGC 2974. The large pitch angle of the gaseous arms ob- 
served in the inner 2^^of NGC 2974 is reminiscent of the 
straight /curved gaseous and dust lanes seen in large-scale 
bars (e.g. Athanassoula 1992). Such structures correspond 
to bar-driven shocks inside the Corotation Resonance of the 
bar (CR hereafter), the dissipative component being driven 
towards the Inner Lindblad Resonance (ILR hereafter) along 
the leading edge of the bar. This obviously requires the pres- 
ence of an Inner Lindblad Resonance (Athanassoula 1992, 
Maciejewski et al. 2002), therefore setting some constraints 
on the central mass profile. 
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Maciejewski et al. (2002) emphasized the characteristic 
features of the gas flow in a single bar. One of the main 
characteristics is the abrupt change of velocity when the gas 
crosses the principal shocks mentioned above, flowing from 
the traiUng to the leading side of the arms. This could qual- 
itatively explain the strong streaming seen in Fig. [13 and 
the large width of the emission line on the trailing side of 
the spiral arms in NGC 2974. Maciejewski et al. (2002) also 
observe gas being funelled towards the 4/1 resonance (or Ul- 
tra Harmonic Resonance, hereafter UHR) along convergence 
regions. This can create gas concentrations at the end of the 
principal shocks, and in the case of NGC 2974 could account 
for the abrupt change in the pitch angle at ^ 2" . This hy- 
pothesis of the presence of an inner bar in NGC 2974 is 
examined in the next section, using the formalism described 
in Emsellem et al. (2001). 



4.4 The bar model 

In Section ^3 we obtained a mass model assuming axisym- 
metry which reproduces the observed stellar kinematics of 
NGC 2974 in the central 18^^ We can obtain a first and rough 
estimate of the pattern speed of a presumed tumbling 
bar by deriving an "azimuthally averaged" resonance dia- 
gram for NGC 2974 using this mass model, which is shown 
in Fig. 1201 The first striking thing is that the — hi/2 ra- 
dial profile predicted by this model (which does not include 
any central dark mass) inevitably implies the presence of 
a single Inner Lindblad Resonance. As emphasized above, 
the ILR is a prerequisite for the building of shock lanes by 
a tumbling bar. The principal shocks should lie inside the 
radius of the Ultra Harmonic 4:1 Resonance, therefore con- 
straining < 1000 (i2i.5~^ km s~^ kpc~^. We also assume 
that the spiral arms are at least partly outside the ILR, thus 
implying that ^ 500 (^21.5"^ km s~^ kpc~^. 

We then add a bar-like component in the model by per- 
turbing the gravitational potential with a simple second har- 
monic function of the form $(r, 0) — — $2(r) cos 2^, where r 
and are the radius and azimuth in the frame co-rotating 
with the bar perturbation, respectively. For $2(r) (m = 2), 
we use the functional form provided by the /c = 3 bar pro- 
vided by Kalnajs (1976; see also Emsellem et al. 2001). Fi- 
nally, gas orbits were calculated using the epicycle approx- 
imation for a potential perturbed by a bar-like component 
(Lindblad & Lindblad 1994; Wada 1994). We included a 
'damping' term as in Emsellem et al. (2001). The free pa- 
rameters of this model are: the damping coefficient A, the 
pattern speed of the bar Qp, the amplitude ^bar (in km^ s~^) 
and radius rbar (in arcsecond) of the perturbation (see a sim- 
ilar but more complete approach in Wada & Koda 2001). For 
each set of parameters, we need to adjust the angle Otar, in 
the equatorial plane of the galaxy, between the bar and the 
line of nodes. 



4-4'^ Comparison with the observed morphology and 
kinematics 

We examined systematically models with values of Qp within 
the constrained range (500 < Qp < 1000 in units of 
km s~^ kpc~^; see above) with steps of 10 km s~^ kpc~^, 
for different values of A, qhar, and rbar- Only values of rbar 
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Figure 20. Resonance diagram for NGC 2974, derived from its 
MGE model. Profiles of 17, Q-k/2, Q^k/2, Q-k/A, are provided. 
The assumed value for the pattern speed of the bar is shown as 
an horizonthal line with Qp = 700 km/s/kpc. 



between 2^4 and 2^7 can reproduce the change of direction 
of the arm at 2^1; the best fit value is rbar = 2'! 6 (or about 
270 pc at 21.5 Mpc). For values of qhar greater than 2 x 10^, 
the pitch angle of the corresponding arms become too high 
with respect to the observed values. As the pattern speed 
increases, the inner side of the spiral arms moves towards 
the centre. As to Qp, we can find reasonably good qualita- 
tive fits to the observed structures if 500 < Qp < 850 (in 
km s~"^ kpc~"^). 

The gas orbits for a model with Qp = 
700 d2i.5~^ km s~^ kpc~\ ^bar = 1 X 10^ and Obar = 55° 
are shown in Fig. 1211 along with the observed morphology 
of the ionized gas. For this model, the ILR, UHR, CR and 
OLR are at radii of lV4, 3^2, 4^9, and 8^5 respectively. The 
projected major-axis of the bar is then at about 35° from 
the major- axis of the galaxy (for i = 60°). The crowding 
of the streamlines correspond to compression regions: the 
model thus predicts the presence of shocks on the inner side 
of the arms due to the convergence of streamlines. This is 
consistent with the high velocity dispersion observed in the 
TIGER maps (Fig. H^ . Note that this model corresponds 
to a perturbation in the potential of less than 2% and a 
maximum torque of 10%: this definitely corresponds to a 
weak bar (see Fig. 1 of Block et al. 2002). 

A detailed comparison between the observed kinematics 
and the one predicted by this model is however beyond the 
scope of the present study, as it would first require a rea- 
sonable fit of the observed line emission distribution (all ob- 
served maps being luminosity weighted). This density wave 
model is only intended to provide a first qualitative attempt 
at reproducing the observed gas distribution and kinemat- 
ics. The existing data are not sufficient to really constrain 
the input parameters of the static model we developed: we 
can only provide rough estimates of the main parameters 
such as: the pattern speed Qp = 700 ± 100 km s~^ kpc~^ 
the CR radius Rcr = 4.9 
and the line of nodes 



-1.4 



Q g, the angle between the bar 
bar = 55° ± 20°. It may be possi- 



ble to narrow the range of possible values for the pattern 
speed of the wave, as its value determines the location of 
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Figure 21. Gas orbits in a model including a bar perturbation in NGC 2974 (bottom panels) as compared to the observed WFPC2 
Hq;+[Nii] image (top panels). The left-hand panels are representations on the sky plane, the right-hand panels correspond to deprojected 
views (assuming an inclination angle of 60°). The crowding of the streamlines reveal a two-arm spiral resembling the observed gaseous 
structure. The line of nodes is horizontal and the major-axis of the bar is Oiyar = 55° away from it (indicated by the dashed line in the 
lower two panels). 



the resonances, but this would require a full hydrodynam- 
ical simulation. A deeper high-resolution Ha+[Nii] image 
would also help to follow the gaseous distribution to greater 
galactocentric distances. 



5 DISCUSSION 

5.1 Dynamical models for NGC 2974 

CvdM94 assumed the presence of a hidden thin stellar disc in 
NGC 2974, contributing a significant 7% to the total visible 
mass of NGC 2974. In the present study however, we solely 
relied on the observed surface brightness (in the / band) 
to build our photometric model, and did not include any 
additional disc component. Our dynamical modeling shows 
that there is no need for the presence of such a hidden com- 
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ponent (see Table |21> to fit the observed stellar kinemat- 
ics. We should however emphasize that the detailed charac- 
teristics of the individual gaussian components provided by 
the MGE formalism are (also) somewhat model-dependent, 
and should be taken with some level of caution. We de- 
rived the corresponding two-integral distribution function 
(strictly speaking, only the even part is constrained by our 
mass model) thus providing a self-consistent model. 

We note that Jeans models do not guarantee that the 
distribution function is positive everywhere. They are gener- 
ally mainly used to constrain the first two velocity moments. 
They should therefore be taken with caution, and should 
only be used, in our view, as a useful tool to provide first 
approximations on global parameters (e.g. mass-to-light ra- 
tio, inclination). Self-consistent models should be favoured 
to derive LOSVDs, as assuming the shape (e.g., Gaussian) 
of the contribution of individual components (bulge, disc) 
may lead to misleading results. 

In our modeling of the stellar kinematics, we made the 
assumptions of axisymmetry and a constant mass-to-light 
ratio, and neglected the effect of dust extinction. The cen- 
tral gaseous kinematics (see Sect. 14.4^ . and the possible stel- 
lar rotation along the minor-axis seem however to invali- 
date the hypothesis of axisymmetry. We also know from the 
i75'T/WFPC2 images that dust is indeed present in the cen- 
tral kpc of NGC 2974, although its distribution is localized. 
The fact that a simple two-integral model fits rather well 
the observed kinematics might just then be due to a lack 
of two-dimensional coverage in the outer part, and to the 
lack of spatial resolution in the central few arcseconds. New 
integral field data acquired with a larger field of view (e.g. 
SAURON on the 4.2-m WHT), better spatial resolution and 
signal- to- noise ratio should thus permit one to further con- 
strain the stellar dynamical structure of NGC 2974. 

5.2 The central dark mass 

If we now interpret the large line width of the central peak 
of the ionized gas (a = 460 km s~^) in terms of a rotat- 
ing Keplerian disc, it implies^ an enclosed mass of about 
5 10^ ^0- This is to be compared with the visible mass 
of 4.7 10'^ ^21. 5 Mq from the MGE model within this ra- 
dius, which would thus be inconsistent with the upper limit 
for a central dark mass of 2.5 10^^21.5 M© fSect.lOl. This 
was however derived under the assumption of a two-integral 
dynamical model, which is not valid if a tumbling bar is 
present. Also we would definitely need data at very high 
spatial resolution to probe the central peak, and check our 
hypothesis of a Keplerian rotating disk. The only instrument 
currently available to help unravel the dynamical nature of 
the innermost region of ionized gas in NGC 2974 is the Space 
Telescope Imaging Spectrograph (STIS) aboard HST. 

5.3 The bar in NGC 2974 

The bar model presented in Sect. 14.41 provides a good repre- 
sentation of the observed gaseous structure and kinematics. 
Within this picture, we can roughly estimate the expected 

^ This depends on the - unknown - radial distribution of the gas 
within the central 0(^1. 



extent of the bar by assuming it is a factor of 1.2 smaller than 
its CR radius: this gives a radius of 4^1 or ^ 425 pc. The 
presence of strong shocks in the gaseous spiral arms then 
implies some dissipation, with a loss of angular momentum 
in the case of the observed trailing spiral in NGC 2974. At 
the scale of the ILR radius (1^4 or 146 pc), the dynamical 
timescale is short, of the order of 6 10^ yr (the rotation pe- 
riod of the suggested bar being 8.7 10^ yr). We therefore 
expect a significant time evolution of the observed structure 
within a few 10^ yr. In this context, the static model based 
on the epicyclic approximation presented in Sect . 14.41 should 
be used with caution, and can only serve as a guideline for 
future full hydrodynamical simulations. 

The gaseous response to a tumbling bar for a potential 
with a single ILR was examined by Fukuda, Wada & Habe 
(1998) via numerical simulations. Their corresponding run 
"Bb" indeed produces trailing spiral shocks resembling the 
ones we observe in NGC 2974. Their simulations suggest 
that fueling is rather inefficient in this case: only gas origi- 
nally within a radius about twice the ILR falls into a central 
structure much smaller than the ILR (roughly 70 pc in the 
case of NGC 2974). The central peak of ionized gas, which 
represents a reservoir of a few 10^ <i2i.5^ (see Sect. l3.2"T)) 
within a radius of less than 0^.^06 (6.25 pc), could be the re- 
sult of such an accretion. 

We finish the discussion on a perhaps more provoca- 
tive note. Considering the relatively low contrast of the pre- 
sumed bar-like perturbation in NGC 2974, and the dust ex- 
tinction, such a structure is extremely difficult to detect us- 
ing visible broad-band photometry, even with the help from 
two-dimensional stellar kinematics. Our result was only pos- 
sible because of the availability of a deep, high-resolution 
image of emission-line gas which showed the spiral structure 
(Figs.ElandEl and the simultaneous use of two-dimensional 
integral-field spectroscopy. This brings up the question of 
how many such bars might be hidden in other early-type 
galaxies? We suggest that these might be more common 
than thought previously, and may provide the possibility 
to funnel gas down to a scale of a few tens of parcsecs, a 
step further towards feeding the AGNs, which are commonly 
found in early- type disk galaxies. In this context, emission 
line gas may certainly be a more helpful tracer than the dust 
features in probing the central structure of nearby galaxies, 
but this then requires both high spatial resolution and two- 
dimensional coverage. 



6 CONCLUSIONS 

In this final Section, we summarize the main results of this 
study. 

• We have presented integral field spectroscopy obtained 
with the TIGER spectrograph which allowed us to probe the 
stellar and gaseous distribution in the central 500 pc of 
NGC 2974. We have applied an iterative scheme to disen- 
tangle the relative contributions of the gaseous emission and 
stellar absorption lines from the obtained datacubes. 

• The stellar line indices do not reveal any significant 
metallicity gradient in the central 500 pc. 

• The ionized gas distribution exhibits a highly con- 
trasted two- arm spiral structure extending up to about 2^^, 
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as well as a central unresolved peak of emission, representing 
a total gas mass of ^ 6.8 10^ M©. 

• The gas kinematics shows strong evidence for non- 
circular motions, with complex and broad line profiles on 
the inner side of the spiral arms. 

• We have detected a central, unresolved line-emitting 
component with a ^ 460 km s~^, both in Ha, and in forbid- 
den lines ([Nil] and [Sii]). Assuming this corresponds to an 
inner Keplerian disc within the central unresolved peak ob- 
served in the narrow-band (Ha+[Nii]) if5'T/WFPC2 image, 
this corresponds to an equivalent mass of about 5 10^ ^0- 

• We built a dynamical model which provides a very good 
fit to the TIGER (and available long-slit) stellar kinemat- 
ics without requiring the addition of an extra component 
such as a disc or a central dark mass. An upper limit of 
2.5 10^^21.5 M0 for the mass of a potential supermassive 
black hole is derived under the assumption of a two-integral 
distribution function: a more general dynamical model could 
however relax this upper limit significantly. 

• A simple density wave (m = 2 bar) model provides 
a good explanation for the observed gas morphology and 
kinematics. Our mass model then inevitably implies the ex- 
istence of a strong Inner Lindblad Resonance, which in turn 
leads to the presence of strong trailing shocks. 

• The best-fit bar model requires a quadrupole pertur- 
bation in the potential of only about 2% (with a maxi- 
mum torque of 10%), with a pattern speed of = 700 =b 
100 km.s~^ kpc~^. This leads to radii for the Corotation 
and Inner Lindblad Resonance of 4^9 and 1^.^4 respectively. 

• Considering the low contrast of the required bar-like 
perturbation, its small size (270 pc in radius), and the pres- 
ence of patchy dust in NGC 2974, such a structure would be 
difficult to detect photometrically in the visible. The pres- 
ence of such inner bars may therefore be more common in 
early- type disk galaxies than generally thought. Still, the 
inner bar in NGC 2974 is of sufficient strength to drive gas 
from a few hundreds of parsecs inwards to the central tens 
of parsecs. 
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